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We present a measurement of the shapes of 6-jets using 300 pb _1 of data obtained with the 
upgraded Collider Detector at Fermilab (CDF II) in pp collisions at center of mass energy = 1.96 
TeV. This measurement covers a wide transverse momentum range, from 52 to 300 GeV/c. Samples 
of heavy-flavor enhanced jets together with inclusive jets are used to extract the average shapes 
of 6-jets. The 6-jets are expected to be broader than inclusive jets. Moreover, 6-jets containing a 
single 6-quark are expected to be narrower than those containing a 66 pair from gluon splitting. 
The measured 6-jet shapes are found to be significantly broader than expected from the pythia and 
herwig Monte Carlo simulations. This effect may arise from an underestimation of the fraction 
of 6-jets originating from gluon splitting in these simulations. The jet shape distributions provided 
in this paper could be compared to any full Monte Carlo simulation and could be used to further 
constrain the various parameters. 
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I. INTRODUCTION 



The measurement of jet shapes allows a study of the 
processes that occur between the initial hard interaction 
and the collimated flow of hadrons observed experimen- 
tally [I]. The internal structure of jets is dictated princi- 
pally by the multiple gluon emissions from the primary 
parton. Multi- gluon emission involves high order QCD 
processes that are complicated to calculate, and vari- 
ous parton shower models are used to implement this 
in Monte Carlo (MC) simulation. In addition to this, a 
good understanding of the hadronization and fragmenta- 
tion processes is needed in order to compare simulation 
results with what is observed in data. For heavy flavor 
jets, the decay of the heavy hadrons must also be cor- 
rectly modeled. Jet shapes, that describe the transverse 
momentum distribution inside jets as a function of dis- 
tance from the jet axis, are defined in detail in Sec. |IV| 
They are interesting distributions to measure in order 
to study the overall decay structure leading to the ob- 
served jets. Moreover, the underlying event, an impor- 
tant component of any hadronic collision, contributes to 
the overall jet shapes. The underlying event comprises 
initial- and final-state radiation, multiple parton interac- 
tions, and beam-beam remnants [2]. 
In this paper, the jet shapes for 6-quark jets are mea- 
sured using the data collected with the upgraded Col- 
lider Detector at Fermilab (CDF II) [3J. The results 
are compared to two different leading order Monte Carlo 
models; the first model is a tuned version of pythia 0], 
the so-called pythia Tune A [2] [5], and the second is 
herwig [6]. 



These models are described in Sec. VI Jet 
shapes are known to be sensitive to whether the initial 
hard-scattered parton is a quark or a gluon; it is also ex- 
pected that jet shapes are sensitive to the initial quark 
flavor. In the case of heavy flavor jets, the shapes are ex- 
pected to be sensitive to the relative contributions of the 
different production mechanisms. The b and the b quarks 
from gluon splitting are expected to be most often inside 
the same final jet [7], leading to significantly broader jet 
shapes than for jets originating from flavor creation. The 
fraction of gluon splitting events is an important param- 
eter for the tuning of Monte Carlo simulations. 
The inclusive jet shapes have been previously measured 
at CDF II [8 . The measured jet shape variables are 
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found to be well described by pythia Tune A MC sim- 
ulation. The data were found to not agree as well with 
the HERWIG MC simulation. 

Charm jet shapes were studied by the HI collaboration 
in photoproduction at the Hadron-Electron Ring Accel- 
erator at DESY (HERA) [9 j . Deviations of pythia from 
measured data were observed in the region where the 
gluon splitting to cc pair is expected, from simulation, to 
contribute significantly. The study found no deviations 
in inclusive di-jet photoproduction. 

Some measurements investigating correlations between 
bb pairs have shown indications of an enhancement over 
the LO expectations of the contribution from gluon split- 
ting. The CDF collaboration [10] measured bb azimuthal 
correlations in pp collisions at 1.8 TeV by requiring two 
jets to be heavy-flavor tagged in semileptonic 6-hadron 
decay and analyzing the azimuthal angle difference 
between these two jets. The transverse momentum range 
was lower than that reported in the present paper. The 
minimum transverse momentum, px [111 , for e or /i lep- 
tons was 8 GeV/c, from which a minimum 6-hadron pt 
of 14 GeV/c was inferred. From the enhanced yield at 
small A(f>, it was concluded that the contribution of gluon 
splitting in heavy flavor had to be roughly doubled over 
expectations from the leading order pythia Monte Carlo 
models. The predictions from a next to leading order 
Monte Carlo simulation show good agreement with the 
data. 

Similar conclusions were reached by the D0 collabora- 
tion |12j in a study of bb production cross section and az- 
imuthal correlation using single muon and dimuon sam- 
ples. The &-hadron pt range was 6 < p\ < 30 GeV/c. 
They used a next to leading order (NLO) Monte Carlo 
simulation, and found good agreement with data in the 
azimuthal correlation between the two muons. The NLO 
calculation enhances the bb production in the A<f> < 1 re- 
gion substantially over the expectations of leading order, 
which has no gluon splitting in the matrix element. 
A third measurement has recently been reported by 
CDF [13 for jet transverse energies Ex > 35 GeV [TTj . 
The heavy flavor content was tagged by secondary ver- 
tices, and A(j> was defined as the azimuthal angle between 
the two secondary vertices. Like the other two studies, 
an enhancement was observed for A(f> < 1, above the 
expectations of pythia and herwig Monte Carlos sim- 
ulations. 

The shape of 6-quark jets provides another method of 
studying heavy flavor production mechanisms. This anal- 
ysis aims to investigate if the fraction of fr-jets originating 
from gluon splitting, as well as its evolution with px, is 
well described in the Monte Carlo models. This is par- 
ticularly important for extrapolations to higher energies, 
such as at the LHC. The comparison of measured vari- 
ables related to the internal structure of 6-jets, such as 
jet shapes, to Monte Carlo simulations is sensitive to the 
global effect of combining models for 6-quark production, 
fragmentation, hadronization, and 6-hadron decay. 
This paper presents a measurement of the integrated b- 
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jet shape variable in inclusive jet production in pp col- 
lisions at y/s — 1.96 TeV. The data were obtained us- 
ing the CDF II detector. The integrated luminosity of 
the data sample, collected from March 2002 to August 
2004, is about 300 pb" 1 . Jets are reconstructed using 
the MidPoint cone algorithm with a cone size of 0.7. Two 
data sets are used for this measurement: a sample of jets 
from inclusive jet production, containing only a negligi- 
ble amount of &-jets and a sample where the heavy flavor 
fraction has been enhanced by requiring a reconstructed 
secondary vertex inside the jet that comes from the decay 
of a heavy hadron. The integrated fr-jet shape variable 
for four different pt regions covering the range from 52 
to 300 GeV/c are extracted from these samples after cor- 
recting for the biases introduced by the tagging as well 
as detector effects, as described in Sec. VII These cor- 
rections are obtained from Monte Carlo samples which 
are passed through a full CDF detector simulation based 
on GEANT3 [H] where the GFLASH [T5] package is used to 
simulate the energy deposition in the calorimeters. The 
results are compared to the leading order Monte Carlo 
predictions from pythia Tune A and herwig. 
This paper is organized as follows. In Sec. [IT] we present 
a brief description of the most important detector sub- 
systems for this analysis. The jet algorithm used to re- 
construct the jets used in this analysis, MidPoint, is de- 
scribed in Sec. IIIII Section ITVl defines the function used 
to describe the jet shapes. Section [V] describes the data 
samples that are used for this analysis along with the 
event selection. A description of the different MC mod- 
els used in this analysis, both in the extraction of cor- 
rection factors and for comparisons of the final results, 
is to be found in Sec. IVII In Sec. IVIII the method used 
to extract the 6-jet shape variable from a sample of in- 
clusive jets and from one of heavy-flavor enhanced jets is 
presented. The systematic uncertainties are presented in 
Sec. |VIII| followed by the results in Sec. |IX| Conclusions 
are summarized in Sec. |Xj 



II. THE CDF II DETECTOR OVERVIEW 

This section presents the CDF II detector, a general 
purpose detector with azimuthal and forward-backward 
symmetry. It is composed of independent subsystems de- 
signed for distinct tasks related to the study of pp inter- 
actions. The two most relevant systems for this analysis 
are the tracking detectors and the calorimeters. A com- 
plete description of the CDF II detector can be found 
elsewhere [3]. 

The tracking system consists of a large open-cell drift 
chamber and silicon microstrip detectors. These lie in- 
side a superconducting solenoid coaxial with the beam 
which produces a magnetic field of 1.4 T. The fiducial 
region of the silicon detector covers the pseudorapidity 
range \rj\ < 2; this subsystem, the closest to the beam- 
pipe and with the finest segmentation, is used for recon- 
structing displaced vertices. The drift chamber measures 



the momentum of charged particles up to \rj\ < 1. Seg- 
mented sampling electromagnetic and hadronic calorime- 
ters, arranged in a projective tower geometry, surround 
the solenoid magnet and measure the energy flow of in- 
teracting particles in the region \rj\ < 3.6. The central 
barrel calorimeters cover the region |?7| < 1, the most 
relevant region for this paper. The segmentation in 4> 
and r\ for both the central electromagnetic and hadronic 
calorimeters is 15° and 0.11, respectively. The transverse 
momentum associated with a given calorimeter tower is 
obtained by assuming that the total tower momentum is 
given by the tower energy; the direction of the momen- 
tum vector is taken to be parallel to the vector linking the 
primary interaction and that tower; the transverse mo- 
mentum is the projection of the momentum vector onto 
the plane perpendicular to the beam direction. Finally, 
a three-level trigger system [16J is used to select events 
online, as described in the section on event selection. 



III. JET RECONSTRUCTION 

In this section, the reconstruction algorithm used to re- 
construct the jets used in this analysis is described. Jets 
used in this analysis are reconstructed using the Mid- 
Point cone algorithm [T7] [TH] , 

Before any jet algorithm is run, the electromagnetic and 
hadronic sections of each calorimeter tower are combined 
into physics towers. Each section is identified with the 
vector joining the primary vertex of the interaction and 
the section's geometrical center. The four- vector mo- 
mentum components, P = (p x ,p y , p z ,E), of each physics 
tower are then computed using the four-momentum sum 
of its electromagnetic and hadronic components; only 
towers with transverse momentum above 0.1 GeV/c are 
considered for jet reconstruction. No contributions for 
the mass of particles are included in this calculation. 
The MidPoint jet algorithm is then run over these physics 
towers. Each physics tower with px > 1 GeV/c is used to 
define a seed around which a jet can be formed. Starting 
from the seed in the event with the highest px, a cone is 
drawn around each seed. The radius of this search cone 
is half the jet cone size, measured in (y,<^)-space |llj . 
The sum of the 4-momenta of all physics towers inside 
the cluster defines the 4-momcntum of the cluster. The 
rapidity, y and <p of the cluster are computed as 



y 



2 y E cl 



Pi 



^tan- 1 ^). 



(1) 



Starting from these clusters, the energy- weighted cen- 
troid including all contributions from the towers within 
the cluster is computed. This new point is then used 
as the center for a new cluster. This procedure is re- 
peated until a stable solution is found. This occurs when 
the geometrical centre of the cluster is aligned with its 
energy-weighted centroid. In the next step, the midpoint 
between each pair of stable cluster centroids, separated 
by less than twice the jet cone radius, is added to the list 
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of cluster centroids. The clustering algorithm is again 
iterated until the new set of stable clusters is found. Fi- 
nally the cluster size is increased from R/2 to the jet cone 
size, R. At this point, if two jets overlap, the momentum 
sharing is considered; if the fraction of the momentum 
of the jet that overlaps with another jet is larger than 
75%, the jets are merged. Else the towers are associated 
to the jet whose center is closest. This jet algorithm can 
be applied in a similar way to final state hadrons in MC 
generated events instead of the physics towers. 
The size of the jet cone is chosen to be 0.7 for this anal- 
ysis in order to be comparable to the inclusive jet shape 
measurement from CDF [5] . A larger jet cone will result 
in a larger fraction of gluon jets where the gluon splits 
into a qq pair and both quarks end up inside the same 
jet. A smaller jet cone might result in this topology be- 
ing reconstructed as two jets, depending on the amount 
of parton showering. 

IV. JET SHAPE VARIABLE DEFINITION 

This section presents the jet shape variable used in 
this analysis. Jet shapes are defined as the distribution, 
as a function of the distance away from the jet axis, 
of the fractional transverse momentum inside the jet; 
jet shapes therefore measure the fraction of the total 
jet transverse momentum inside a given radius in 
(y,</>)-space from the jet axis. The integrated jet shape 
variable, ^(r/R), is defined as the fraction of total pr 
inside the jet cone of radius R that is carried by the 
particles in a sub-cone of radius r. This quantity can 
be computed at the hadron level, using all final state 
hadrons, in the Monte Carlo simulations and by using 
calorimeter energy deposits or charged particle tracks 
for both the data and the full simulation. The measured 
quantity is the average integrated jet shape, which is 
computed over an ensemble of jets. This quantity is 
expressed as 

where pt(0 — > r) is the scalar sum of the transverse mo- 
menta of all objects inside a sub-cone of radius r around 
the jet axis. The integrated shapes are by definition 
normalized such that *ff(r/R — 1) = 1. By definition 
^(0) = 0. In this analysis, we do not consider particles 
outside the jet cone radius. In Monte Carlo simulation, 
6-jets are defined as jets which have at least one 6-quark 
inside the jet cone. 

V. EVENT SELECTION 

In this section the online event selection criteria (trig- 
gers) used for this measurement are described. This is 
followed by a description of the method used to enhance 



the heavy-flavor fraction of jets by requiring there to be 
a reconstructed displaced vertex inside the jet cone. The 
selections (cuts) applied to the events are then introduced 
and finally the correction factors applied to the jet pt are 
described. 

This paper presents results for central jets, \y\ < 0.7, in a 
Pt range from 52 to 300 GeV/c. Four different datasets 
are used. Events were collected that satisfy the condi- 
tions required by the inclusive jet trigger, with different 
minimum transverse energy thresholds for the different 
datasets. Each dataset is defined by a unique trigger 
path that has unique requirements at each of the three 
stages of the trigger (Trigger Levels LI, L2 and L3) . 
The inclusive jet triggers select events based solely on 
the energy deposits in the calorimeters, with four differ- 
ent thresholds on the jet transverse energy, Et (see Ta- 
ble |l|. Due to limited bandwidth, only a fraction of the 
events in each trigger path is accepted [19] ■ In the first- 
level trigger, a single trigger tower [5U] with Et above 5 
or 10 GeV, depending on the trigger path, is required. 
In the second-level trigger, calorimeter tower clusters are 
formed around the selected trigger towers. Starting from 
the initial first-level tower, adjacent towers are associated 
to it if their Et is above 1 GeV, thus forming clusters 
of calorimeter towers. The events are required to have 
at least one second-level trigger cluster with Et above 
a given threshold, which varies between 15 and 90 GeV 
for the different trigger paths. In the third-level trigger, 
jets are reconstructed using the CDF Run I cone algo- 
rithm [21] , and the events are required to have at least 
one jet with Et above selected thresholds between 20 
and 100 GeV. 

The inclusive jet triggers are not fully efficient at the 
trigger threshold with respect to offline reconstructed 
jets. To avoid any trigger bias, events in each dataset 
are only considered when the trigger efficiency is above 
99%. These trigger efficiency thresholds are shown in 
the last column of Table [I] where the jet transverse mo- 
mentum is quoted after corrections for detector effects, 
as discussed later in this section. The lowest pt for this 
analysis, 52 GeV, is determined by the trigger threshold 
for the lowest pt inclusive jet trigger used. 
The inclusive jet datasets are dominated by light quark 



Trigger Path 


LI 


L2 L3 


99% trig. eff. 




Tower 


Cluster Jet 


Offline 




E T [GeV] 


E T [GeV] E T [GeV] 


Pt [GeV/c] 


JET 20 


5 


15 20 


52 


JET 50 


5 


40 50 


80 


JET 70 


10 


60 70 


104 


JET 100 


10 


90 100 


142 



TABLE I: Summary of trigger paths and trigger thresholds 
used in each of the three CDF trigger levels. The last column 
shows the final offline cuts applied to the corrected jets after 
the 99% trigger efficiency requirement has been applied. 



and gluon jets. The fr-jet content is enhanced by using 
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a secondary vertex tagger [55]. This tagging algorithm 
exploits the long lifetime of the 6-hadrons. Due to the 
large relativistic boost of the particles, most 6-hadrons 
travel a few millimeters before decaying. The tagging al- 
gorithm is based on the reconstruction of a displaced, or 
secondary, vertex using the reconstructed charged parti- 
cle trajectories, tracks, contained within a cone of 0.4, 
in (y, 0)-space, around the jet axis. Despite the jet cone 
size being 0.7, the cone for finding displaced tracks, the 
tagging cone, is kept at 0.4 because the direction of the 
heavy-flavor hadron and its decay products tend to be 
close to the direction of the jet axis. Using a larger cone 
size for the tracks would lead to an increase in misiden- 
tified secondary vertices. The tagging algorithm imple- 
mentation is relatively complex and is described in |22j ; a 
brief overview is given here. Before attempting to recon- 
struct a displaced vertex, several quality cuts are applied 
to the tracks inside the tagging cone in order to reduce 
the number of tracks that are not well measured. Tracks 
are then ordered according to quality criteria, including 
their distance of closest approach to the primary vertex, 
d,Q. Starting from the track with the highest quality, an 
attempt at creating a displaced vertex is made with the 
next best track. The attempt succeeds or fails depend- 
ing on the quality of the fit. If the attempt succeeds, 
all other tracks whose do significance [23J with respect 
to the displaced vertex is below a predefined limit are 
attached to it. If the attempt fails, the highest quality 
track is associated to the third best track and the process 
starts over. If no vertex is found with the highest qual- 
ity track, the second highest quality track is considered 
and attempts are made to reconstruct a displaced ver- 
tex with the next best track. This loop continues until 
either a good displaced vertex is found or there are no 
more pairs of tracks that pass the selection criteria. If a 
good displaced vertex is found, a final cut is applied on 
the significance of the two-dimensional projection along 
the jet axis of the distance between the primary and sec- 
ondary vertex locations [21] . The main limitation of this 
tagging method is the fact that c-hadrons can have decay 
lengths similar to those of 6-hadrons. Another limitation 
is the finite resolution of both the primary and secondary 
vertex location. 

The following list summarizes all the cuts applied to the 
events in the samples: 

One, and only one, primary vertex with z-component, 
z vtx , in the region |z v tx| < 50 cm to ensure good sec- 
ondary vertex reconstruction and to reject events with 
multiple pp interactions; missing transverse energy sig- 
nificance of the event, Pt/\/Et [25], must be smaller 
than a given threshold, varying from 3.5 to 7.0 GeV 1 / 2 
depending on the trigger used, in order to reject cosmic 
rays that enter the detector as well as to reject beam in- 
duced backgrounds; 

1 2/jet | < 0.7 where the secondary vertex tagging algorithm 
is best understood; 

for the 6-jet enriched samples only, jets are required to 
have a secondary vertex tag [22] , 



An average jet correction is applied to correct the jet 
Pt to the hadron level, i.e. to remove any bias due to 
the detector. It is calculated by matching, in the Monte 
Carlo simulation, hadron and calorimeter level jets in 
(y ~ </>)-space. The jet pt is obtained by running the 
jet reconstruction algorithm over the final state hadrons 
and energy deposits in the calorimeter towers, respec- 
tively. The hadron level jet px is plotted as a function of 
the calorimeter jet pr- The scatter plot obtained is fit- 
ted, using a chi-squared minimization, to a fourth order 
polynomial which is then used to correct the transverse 
momentum of each measured jet in data. The increase in 
Pt due to this correction is on the order of 20% at low pr 
and 10% at higher pt- The binning used for this analysis 
and for all plots shown is in corrected jet pt, referred to 
in the remainder of this paper simply as the jet px- This 
is an absolute correction to the jet that does not affect 
the shape of a given jet in the sample, only its total pt 
and hence which px bin it belongs to. The corrections to 



the jet shape variables are discussed in Sec. VII E 



VI. MONTE CARLO MODELS 

This section introduces the two leading order MC mod- 
els to which the measured 6-jet shapes are compared: 
PYTHIA Tune A and herwig. In both MC models, 
the parton shower is implemented to leading log. For 
both MC models, inclusive dijet samples were produced 
(mscl = 1) using the CTEQ5L parton distribution func- 
tions [25] , 

PYTHIA Tune A is based on PYTHIA v6.203 g] which has 
been tuned to the CDF Run I underlying event. It is 
worth noting that the underlying event tuning is found 
to be important for a good description of the inclusive jet 
shape [8j. This underlying event tune is widely used in 
CDF analyses [8]. The final-state parton shower model 
of pythia used for this analysis is carried out in the 
jetset part of the code [27]; this time-like evolution 
is computed using splitting functions and an ordering in 
virtuality scale, Q 2 . The splitting functions, P a ^bc{ z ), 
describe the probability that a parton a splits into two 
partons b and c where parton b carries away a momen- 
tum fraction z of the initial parton. The virtuality scale 
in pythia is equal to the four-momentum squared of the 
branching parton. The ordering in virtuality scale means 
that for initial state radiation, Q 2 increases as the pro- 
cess unfolds from the incoming primary parton as it ap- 
proaches the hard collision. For final state radiation, Q 2 
decreases going from the energetic outgoing parton to the 
final shower. Decreasing Q 2 in the final state is similar 
to the decreasing angular ordering used by herwig. The 
parton shower is terminated at Q of 1 GeV/c for QCD 
branchings. For the non-perturbative fragmentation and 
hadronization processes, our version of pythia uses the 
default Lund string model, as implemented in jetset. In 
this model, the long-range confinement forces are allowed 
to distribute the energies and flavors of a parton configu- 
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ration among a collection of primary hadrons which sub- 
sequently may decay further. 

The modeling of the decay of 6- and c-hadrons is based 
on the program QQ [55], the decay tables of which are 
periodically updated with the latest experimental infor- 
mation. When pythia produces a 6-hadron in a jet, QQ 
is called to handle the decay. For some decay modes the 
exclusive final states are only obtained after fragmenta- 
tion of the quarks, which requires a return to pythia. 
The second model to which the experimental results are 
compared is herwig v6.505 [BJ. In the perturbative par- 
ton showering process, the main difference with respect 
to pythia is the use of angular ordering of successive 
emissions, which simulates the color flow of the subpro- 
cesses. For final state radiation, the angle of the radiated 
gluon relative to the parent parton direction decreases as 
the process unfolds. The maximum angle is determined 
by the elementary subprocess and is due to interference 
among gluons. herwig uses the cluster model for the 
fragmentation into hadrons [29 . This model, which is 
independent of the initial hard process and of the energy, 
is intended to disrupt as little as possible the event struc- 
ture established by the parton showering. Color-neutral 
clusters are formed that decay into the observed hadrons. 
The decay of heavy flavored hadrons is treated in HER- 
WIG in the same way as for all other unstable particles, 
according to the exponential decay law using the mean 
lifetime of the particle. This version of herwig does not 
include multiple parton interactions, which is thought to 
be the main reason why the inclusive jet shapes agree 
better with pythia Tune A than herwig [51. 



VII. B-JET SHAPE UNFOLDING 

In this section, the analysis methodology is presented 
which makes it possible to extract the shapes of 6-jets 
from samples of inclusive and tagged jets. First the un- 
folding equations are described, then in Sec. |VII A| the 
effect on these equations of the fraction of 6-jets that 
contain more than one 6-quark is discussed. The next 
sub-sections present how the different components of the 



unfolding equations are measured; Section fyTI B| presents 
the measurement of the inclusive and tagged samples; 
Section |VII C| presents the method used to extract the 
6-jet fraction in the tagged samples; the method used to 
account for the bias to the jet shapes due to tagging is 
presented in Sec. |VIID| finally the corrections applied to 
the measured 6-jet shapes in order to obtain a detector 
independent result are presented in Sec. |VII E| 
The jet shape variables, defined in equation [2] are 
computed in data using the calorimeter towers with 
Pt > 0.1 GeV/c. The use of calorimeter towers with 
Pt > 0.5 GeV/c as well as the use of tracks with 
Pt > 0.5 GeV/c are considered in the study of the sys- 
tematic uncertainties. The sample of tagged jets used for 
this analysis does not contain only 6-jets but also back- 
ground jets, called non6-jets, that do not contain any 



6-quarks. The non6-jets contain c-jets, light flavor-jets 
and gluon-jets, where the gluon does not fragment into 
a 66-pair. The purity, p^, is defined as the fraction of 
tagged jets that are 6-jets. The detector level jet shapes 
for tagged jets, ^Jjet> wm * nus ^ e gi yen by 



Pb*£f(r/iZ) 



(l-PbK a r° n V/i?), 



(3) 



where "f^f b ano - ^dot" 01 * are ^ ne 6- & n d nonfo-jet contri- 
butions to the measured tagged jet shape. 
The use of the secondary vertex tagging method to en- 
hance the heavy flavor content of the sample biases the 
measured jet shapes due to the fact that the secondary 
vertex reconstruction algorithm requires jets with clean 
and well defined tracks. A bias term, dependent on the 
distance from the jet axis, r, and on the px of the jet must 
thus be added to correct for this effect without affecting 
the overall normalization of the integrated jet shapes. 
This bias is different for 6-jets and non6-jets. In the case 
of non6-jets, the bias terms also take into account the 
enhanced fraction of c-jets in the tagged jet sample. The 
bias terms, b\>(r/R) and b non b(r/R) for b- and non6-jets 
respectively, are defined in MC simulation such that 



^ h (r/R) = b h (r/R)V h dct (r/R) and 
(r/R) = b nonh (r/R)^ h (r/R), 



- det 
.T.tag nonb 
*det 



(4) 
(5) 



where the \& d e %(r / R) terms represent the detector level 
shapes, as obtained from MC simulation, before any tag- 
ging requirements and ^J tag (r / R) are those obtained from 
MC simulation after the tagging requirement is applied. 
The bias terms are computed separately for each px bin. 
Combining Eq. ^ with the definition of the bias terms 
from Eq. the measured detector level jet shapes for 
tagged jets can be re- written as 



*ll(r/R) - Ph b h (r/R)* h dct (r/R) 



+(l-P b )bnon b (r/RmT(r/R), (6) 

Rearranging Eq. ([6|, the detector level measurement of 
the 6-jet shape is defined as 



dot 



(r/R) 



SyVg) - (1 - Pb )b nonh (r/R)^™ b (r/R) 
p h b h (r/R) 



(7) 

It is also necessary to correct the 6-jet shapes for de- 
tector effects, i.e. to remove all influence of the tracker 
or calorimeters on the measurement. The corresponding 
hadron level correction factors, C had (r/R), are defined 
as 



*Ld(r/fi) = C h * d (r/R)* h det (r/R), 



(8) 



where in MC simulation the shapes are computed both 
at the detector level, ^ det (r/R), and using the final 
state particles, ^f^ ad (r/R). These definitions of the tag- 
ging bias and hadron level correction factors ensure the 
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proper normalization of the jet shapes at each stage of 
the unfolding procedure. They are defined such that 
C had (r/R = 1) = 1, b h (r/R = 1) = 1 and b nonh (r/R = 
1) = 1. 



The final equation used to obtain the hadron level 6-jet 
shape is obtained by combining Eq. ^ and Eq. ^8l. It 
can be written as 



*had(^) 



C had (r/J ^d a >/^) - (l-p h )b nonh (r/R)^ h (r/R) 



p h b b (r/R) 



(9) 



r 



where ^^(t/R) is the measured jet shape for the tagged 
jet sample and *§°? (r/i£) is the measured inclusive jet 
shape, described in Sec. |VIIB| The other parameters of 
this equation are discussed in the following sub-sections. 



We evaluate the systematic effect of this particular choice 
of correction to the fib fraction in Sec. |VIII| 



B. Detector Level Jet Shapes 



A. Single 6-quark jet content 

Many of the distributions used for the extraction of 
the b-jet shapes described in the previous section are 
expected to be different depending on whether the jets 
contain one or two 6-quarks. This section describes how 
this is dealt with in this measurement. The parameters 
used in the unfolding are sensitive to the fraction of 6-jets 
that contain a single 6-quark, /ib- In leading order (LO) 
Monte Carlo, gluon splitting to 66 pairs occurs as part 
of the fragmentation process and not in the matrix cle- 
ment. For most jets where gluons split to a bb pair, both 
6-quarks end up inside the same jet cone [7] . A compar- 
ison between the fraction of 6-jets with more than one 
6-quark inside the same jet cone, 1 — /ib, predicted by 
PYTHIA Tune A and the next to leading order (NLO) 
calculation is shown in Fig. [l] [30] [31]. The maximum 
deviation between the PYTHIA Tune A and the NLO pre- 
diction with the factorization and renormalisation scales, 
/x, equal to ^ is on the order of 0.2. This particular 
choice of /i is motivated by the measurement of the in- 
clusive b-jet cross section [31]. Before calculating any of 
the unfolding factors, the Monte Carlo samples are re- 
weighted to decrease the /ib fraction by 0.2, in order to 
account for this underestimation of the gluon splitting 
fraction. This value was chosen for all pt bins as it cor- 
responds to roughly the expected shift. For each of the 
6-jet shapes obtained from MC simulation: the unbiased 
6-jet shapes, the tagged 6-jet shapes and the hadron level 
6-jet shapes, the shapes are extracted separately for sin- 
gle, \P lb , and double 6-quark jets, ^ 2b . They are then 
combined using the new fraction of single 6-quark jets as 
follows 



* b = (/ib - 0.2)* lb + (1 - (At - 0.2))* 



2b 



(10) 



This new lower /ib fraction is used for the corrections 
to data in the remainder of the analysis, to obtain the 
secondary vertex mass templates, as well as the tagging 
biases and hadron level corrections to the 6-jet shapes. 



In the jet samples considered for this analysis, the av- 
erage multiplicity of calorimeter towers inside the jet in- 
creases slightly as the jet px increases from 16 to 23 tow- 
ers and is on average a little higher for tagged jets than 
for inclusive jets, which range from 13 to 18 towers. The 
multiplicity of tracks that pass all the selection cuts ex- 
hibits the same behavior. The variation is from 16 to 23 
tracks for tagged jets, and 13 to 17 tracks for inclusive 
jets. Jets with a higher pt have more tracks due to the 
fact that the final number of particles is determined by 
the number of times the initial parton radiated another 
parton, before it hadronizes. Jets with higher p? have on 
average a higher initial Q 2 which means a larger range 
over which to radiate more particles. 
The integrated jet shape is computed for each jet i in the 
sample as 



Efe Vt\ x < i) 



(11) 



where the sum is over all towers that satisfy the con- 
ditions of belonging to the jet and that the fractional 
distance to the jet axis, x, must be less than r/R in the 
numerator and 1 in the denominator. The use of this 
denominator ensures the correct shape normalisation for 
each jet. For each p? bin, the average jet shape is com- 
puted from the jet shapes of all jets in that sample 



*(r/i?) =< Vi(r/R) > 



(12) 



The tagged jet shapes, * tag (r/i?), are defined as the av- 
erage jet shapes, measured at the detector level, of all 
tagged jets in the samples. 

Given the very low fraction of 6-jets in inclusive jet pro- 
duction, estimated from the MC to be less than 4%, it is 
possible to approximate the non6-jet shapes to those of 
the inclusive jet shapes, before any tagging requirements. 
The assumption is that ^°? h (r/R) m * d n c c t '(r/i?). The 
difference between these shapes, in pythia Tune A 
Monte Carlo, is negligible with a maximum difference 



10 



of less than 0.5%. No systematic uncertainty is therefore 
related to the use of this approximation. 
Both the inclusive jet shapes, ^f^(r/R) and the tagged 
jet shapes, ^^(r/R), are measured from data in each 
of the four px bins. These detector level jet shapes are 
shown in Fig. [2] for all four px bins. Comparing these 
inclusive and tagged jet shapes shows that there are sig- 
nificant differences in the measured shapes, yielding con- 
fidence that this measurement is possible. But at this 
stage, it is not possible to conclude anything about the 
shapes of 6-jets as the displaced vertex requirement biases 
the observed tagged shapes, as described in Sec. |VIID| 



C. Purity 

The fractions of 6-jets in the tagged jet samples are 
extracted from fits to invariant mass distributions for b- 
and non5-jets, calculated in pythia, from the charged 
tracks forming the secondary vertex. It is not possible 
to reconstruct the full hadron invariant mass mainly be- 
cause of the presence of neutral particles in the 6-hadron 
decays that are not detected in the tracking detectors. 
Nevertheless, the distribution of the invariant mass of the 
tracks associated with the secondary vertex, referred to 
as the secondary vertex mass, is significantly different for 
real heavy flavored jets than for mis-reconstructed light 
flavored or gluon jets. Using the Monte Carlo samples, 
distributions of the secondary vertex masses for tagged 
jets are obtained for each px bin, separately for b- and 
nonfr-jets. The non6-jet distribution is a combination of 
real displaced vertices from c-jets as well as light fla- 
vor jets where a secondary vertex was mistakenly iden- 
tified by the secondary vertex algorithm. As an exam- 
ple, the secondary vertex mass distributions, as obtained 
from PYTHIA Tune A, for the second px bin, from 80 to 
104 GeV/c, are shown in the top two plots of Fig. [3] The 
measured distribution in the data is fitted to the b- and 
nonfr-templates, using a binned y 2 minimization method, 
to find the most probable fraction of jets that are 6-jets, 
as shown in the second plot of Fig. [3] for the second px 
bin. The fit describes the data very well in all px bins. 
Figure [3] (bottom) shows the extracted purity, pb, as a 
function of the pt of the jets as obtained by fitting the 
data with templates from pythia Tunc A on the one 
hand and from HERWIG on the other. 



D. Biases Due to Secondary Vertex Tagging 

The requirement that the jets be tagged by the sec- 
ondary vertex algorithm introduces a bias in the mea- 
sured jet shapes. These biases are different for each pt 
bin, each bin in r and different for b- and non5-jets. The 
bias terms are defined as the ratios, as obtained from 
the Monte Carlo samples, between the tagged and the 



unbiased jet shapes for b- and non6-jets separately 

bhir/R) = ^c(W (13) 

^ b (^) = ||^|. (14) 

The bias term for nonfr-jets takes into account the in- 
creased fraction of c-jets in the tagged jet sample. The 
maximum bias for 6-jets is on the order of 8% and for 
nonfe-jets is on the order of 18%. 

In this case the 6-tagging efficiency is not relevant as we 
are interested in the distortions to the average jet shape 
arising from the tagging and not in absolute cross sec- 
tions. 



E. Hadron Level Corrections 

The hadron level correction factors that remove the 
influence of the detector on the measured 6-jet shapes, 
C had (r/i?), are evaluated from the Monte Carlo samples 
for each bin in r and each bin in px and are defined as 

< n>had/--/p\ _ ^hadMc( r /-ft) /i r\ 

C ir/R) -^uc(r/Ry (15) 

where ^detMcC 7 "/-^) are the Monte Carlo simulated 6-jet 
shapes computed at the detector level, and v I'h ac jMc( r '/-^) 
are the Monte Carlo simulated 6-jet shapes computed 
using final state hadrons. These correction factors are 
on the order of 3% at most. 

VIII. SYSTEMATIC UNCERTAINTIES 

The different sources of systematic uncertainties for 
this measurement are described in this section. 
To account for the sensitivity of the unfolding method to 
the variation of the fxb fraction, the fraction is decreased 
by 0.5. The difference in the measured 6-jet shapes when 
using the fib — 0.5 instead of the default /ib — 0.2, dis- 
cussed in Sec. |VII A] is taken as a systematic uncertainty. 
The samples of tagged jets that are from non6-jets con- 
tain a significant fraction of c-jets; this comes from the 
fact that c-jets have real displaced vertices that are sim- 
ilar to those of 6-jets and often get tagged by the sec- 
ondary vertex algorithm. This is unlike the light jets 
where the reconstructed secondary vertex is not a true 
displaced vertex. The fraction of nonfc-jets in the tagged 
jet sample that are c-jets is taken from the Monte Carlo 
predictions and is found to be between 25% and 50%, 
depending on the px of the sample. From secondary ver- 
tex mass fits separating b-, c- and light-jets into three 
independent templates, a systematic uncertainty of 5% 
is assigned to the c-jet content of the nonfr-jets. This un- 
certainty affects both the secondary vertex mass fit and 
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the tagging bias for non6-jets. 

The c-jets in the non6-jet sample are expected to have 
a somewhat similar behaviour to the 6-jets, in particu- 
lar regarding gluon splitting to cc pairs. Similarly to the 
/ib fraction for 6-jets, one can define the fraction of c-jets 
containing only one c-quark, f\ c . This fraction might also 
be overestimated in the LO Monte Carlo as one would ex- 
pect /ic to increase if /ib were to increase. In order to 
investigate the effect of a possible underestimation of this 
effect, the /i c fraction is decreased by 0.2. This change 
affects both the tagging bias for nonfo-jets as well as the 
templates used in the secondary vertex mass fit. The 
difference in the measured 6-jet shapes when varying the 
/ic fraction is taken as a systematic uncertainty. 
In order to evaluate the effect of using a particular set 
of Monte Carlo models for the fragmentation, hadroni- 
sation and underlying event, the whole unfolding is per- 
formed using herwig Monte Carlo samples instead of the 
PYTHIA Tune A samples. The difference in the measured 
&-jet shapes obtained using these two Monte Carlo sam- 
ples is taken as a systematic uncertainty. 
To gauge the effect on the unfolding procedures of any 
possible mis-modeling in the simulation of the detector 
response, the jet shapes are measured using the charged 
tracks inside the jet cone, instead of the calorimeter tow- 
ers. The measured track jet shapes are unfolded back 
to hadron level using new correction functions obtained 
from MC simulation. The jet direction and transverse 
momentum remain unchanged with respect to the default 
scenario. All tracks with p? > 0.5 GeV/c that originate 
within a cone of 0.7 from the jet axis are considered for 
the measurement of these jet shapes. At detector level, 
the jet shapes defined by tracks tend to be narrower than 
those defined by calorimeter towers. The hadron level 
corrections will thus be different for the measurement 
using tracks than for the one using calorimeter towers 
The difference between the final hadron level jet shapes 
computed using this method and the default one is taken 
as a systematic uncertainty. 

To investigate the accuracy of the simulation of the 
calorimeter response to low px particles, the analysis 
is performed using only calorimeter towers with px > 
0.5 GeV/c and the difference with respect to the nom- 
inal measurement (which uses calorimeter towers with 
Pt > 0.1 GeV/c) is taken as a systematic uncertainty 
and is found to be negligible. 

A 3% systematic uncertainty on the jet energy correc- 
tions is considered that combines the 3% systematic un- 
certainty for inclusive jets |32j with the uncertainty on 
the 6-jet fragmentation that is 0.6% . A variation of 
±15% on the Missing Ex significance is applied. The cut 
on the location of the primary vertex is varied by ±5 cm 
around the nominal cut at 50 cm. These variations are 
all found to have only small effects on the final measure- 
ment. The dependence on the Monte Carlo modeling of 
the secondary vertex parameters was also investigated 
and found to be negligible. 

The total, statistical and systematic uncertainties are 



shown in Fig. [4] for each px bin and r bin. Also shown 
are the various contributions from the dominant effects. 
The dominant sources of systematic uncertainties vary as 
a function of the px bin. These are: 
the difference in the 6-jets shapes reconstructed using 
pythia Tune A and herwig ; 

the difference in the &-jet shapes reconstructed from 
tracks instead of calorimeter towers, i.e. due to the de- 
tector simulation; 
the /ib variation from -0.2 to -0.5; 
the /i c variation by -0.2; 
the jet energy scale. 



IX. RESULTS 

The final results are presented and discussed in this 
section. The measured integrated fr-jet shapes are shown 
in Figs. [5] and [6] as black open squares. The statistical and 
total uncertainties on the measurements are shown; the 
statistical uncertainties are smaller than the points. The 
uncertainty on the measured jet shape is at r/R — 1 
because the integrated jet shape is defined to be exactly 
equal to 1 at this point. The first three bins in r/R thus 
contain most of the information about the broadness of 
the jet shapes. 

The results are compared to pythia Tune A and herwig 
predictions using the default /ib fractions. Figure [5] also 
shows the pythia Tune A predictions for the inclusive 
jet shapes. Reference [5] shows good agreement between 
data and the MC predictions for inclusive jet shapes. 
This plot shows that despite relatively large systematic 
uncertainties, the measurements differ from the inclusive 
jet shape predictions, thus indicating that the jet shapes 
are sensitive to the presence of heavy flavor particles. Jets 
containing 6-quarks appear to be broader than inclusive 
jets. No reasonable change in f-lb could bring the data 
into agreement with the inclusive jet shapes. This plot 
also shows that the agreement between data and LO MC 
simulation is much better using a smaller fraction of jets 
that contain a single 6-quark than the default value. 
The pythia Tune A and HERWIG predictions, when the 
expected distributions if /ib is decreased by 0.2, along 
with the predictions for single and double 6-quark jets, 
are reported in Fig. [6] Single 5-quark jets are predicted 
to be narrower than inclusive 6-jets; double 6-quark jets 
are predicted to be broader. The measured average b- 
jet shapes are between these two curves, narrower than 
double fo-quark jets but broader than single fo-quark jets. 
It appears from these plots that a larger fraction of dou- 
ble 6-quark jets than predicted by the LO MC simulation 
agrees better with the data. The MC method predictions 
for pythia Tune A and HERWIG are shown to be similar 
for all cases considered. 

Figure [TJshows the ratio of the hadron level integrated jet 
shapes from Monte Carlo predictions over the measured 
values for each of the four jet px bins. The light grey 
bands show the total uncertainty on the measurement, 
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whereas the dark grey band shows the statistical uncer- 
tainty. From these plots it is clear that the agreement 
between data and MC simulation is improved in each pt 
bin by decreasing the single 6-quark fraction by 0.2. 
Table |TT| reports the measured 6-jet shapes in each of the 
four pt bins for each of the bins in r. The central values 
are shown along with the statistical and systematic un- 
certainties. 

Another way of looking at these results is to plot the 



52 < p T < 80 GeV/c 
r/R <S> b {r/R)±a st!it ±V s 



0.1/0.7 
0.2/0.7 
0.3/0.7 
0.4/0.7 
0.5/0.7 
0.6/0.7 
0.7/0.7 = 1.00 



0.14 
0.28 
0.42 
0.57 
0.71 
0.86 



0.283 ± 0.010 ± 
0.553 ± 0.010 ± 
0.717 ± 0.007 ± 
0.825 ± 0.005 ± 
0.901 ± 0.003 ± 
0.953 ± 0.002 ± 
1.000 ± 0.000 ± 



0.105 
0.076 
0.068 
0.037 
0.015 
0.006 
0.000 



80 < pt < 104 GeV/c 



r/R 

0.1/0.7 « 

0.2/0.7 « 

0.3/0.7 « 

0.4/0.7 « 

0.5/0.7 « 

0.6/0.7 « 

0.7/0.7 = 



0.14 
0.28 
0.42 
0.57 
0.71 
0.86 
1.00 



tFjr/R) ± ffstgt 
0.336 ± 0.007 ± 
0.565 ± 0.006 ± 
0.710 ± 0.004 ± 
0.817 ± 0.003 ± 
0.898 ± 0.002 ± 
0.957 ± 0.001 ± 
1.000 ± 0.000 ± 



i ^"sys 

0.059 
0.051 
0.039 
0.024 
0.017 
0.006 
0.000 



104 <p T < 142 GeV/c 



r/R 
0.1/0.7 « 
0.2/0.7 « 
0.3/0.7 « 
0.4/0.7 « 
0.5/0.7 « 
0.6/0.7 « 
0.7/0.7 = 



* b (r/E)±a stat ±a s 



0.14 
0.28 
0.42 
0.57 
0.71 
0.86 
1.00 



0.403 ± 0.008 ± 
0.623 ± 0.007 ± 
0.747 ± 0.005 ± 
0.837 ± 0.003 ± 
0.906 ± 0.002 ± 
0.963 ± 0.001 ± 
1.000 ± 0.000 ± 



0.064 
0.024 
0.026 
0.019 
0.010 
0.004 
0.000 



r/R 
0.1/0.7 
0.2/0.7 
0.3/0.7 
0.4/0.7 
0.5/0.7 
0.6/0.7 
0.7/0.7 



142 <p T < 300 GeV/c 



0.14 
0.28 
0.42 
0.57 
0.71 
0.86 
1.00 



<S> b (r/R)±a s ^ 
0.413 ± 0.008 ± 
0.637 ± 0.006 ± 
0.760 ± 0.005 ± 
0.849 ± 0.003 ± 
0.919 ± 0.002 ± 
0.966 ± 0.001 ± 
1.000 ± 0.000 ± 



i & S y S 

0.048 
0.037 
0.020 
0.013 
0.007 
0.008 
0.000 



TABLE II: Integrated jet shapes for 6-jets. The central value 
along with the statistical and systematic uncertainty for each 
Pt and r bin are shown. 



fault /ib fractions as well as the expected distributions 
if /ib is decreased by 0.2. The PYTHIA Tune A predic- 
tions for the inclusive jet shapes as well as the previously 
published inclusive jet shape results are shown in the top 
plot. The rapidity region considered for the inclusive jet 
shape measurement does not include the center-most ra- 
pidity region (\y\ < 0.1). The inclusion or not of this 
region was found in this analysis not to change the value 
of the predictions or the measured values. This figure 
indicates that the evolution of the jet shape with px ap- 
pears to be somewhat flatter for 6-jets than for inclusive 
jets and confirms that the measured b-jet shapes are dif- 
ferent from those measured for inclusive jets. 
The PYTHIA Tunc A and HERWIG predictions for single 
and double 6-quark jets are included on the bottom plot. 
This figure shows that as the px of the jet increases, 
the difference between the pythia Tune A and herwig 
predictions increases, herwig predicts a slightly flatter 
evolution with pt than pythia Tune A. 



X. CONCLUSIONS 

We have reported on a measurement of the b-jet 
shape at the Tevatron collider. Despite the consider- 
able uncertainties of Monte Carlo simulations of this non- 
perturbative process, we show convincing evidence that 
b-jets are broader than inclusive ones. This confirms that 
the jet shape is sensitive to the heavy flavor content. 
The measured 6-jet shapes are significantly broader than 
expected from both pythia Tune A and herwig. One 
possible interpretation is that this effect is coming from 
an underestimation in LO MC simulation of the fraction 
of 6-jets originating from gluon splitting. NLO calcu- 
lations predict a significantly higher rate of 6-jets that 
contain more than one 6-quark inside the jet cone than 
the LO Monte Carlos. Decreasing the relative fraction of 
single 6-quark jets, i.e. increasing the double 6-quark jet 
fraction, by 0.2 in pythia Tune A and herwig leads to 
a better description of the measured 6-jet shapes. This 
decrease is qualitatively consistent with the NLO predic- 
tions with a small factorisation and renormalisation scale, 



^rf, where fi 



These findings are con- 



sistent with a number of other analyses that investigated 
the azimuthal correlations in 66 production. 
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FIG. 1: Fraction of 6-jets that contain more than one 6-quark inside the same jet cone. The pythia Tune A MC simulation 
predictions are compared to NLO calculations for two different hadronisation and factorisation scales. The NLO calculations 
are shown binned with the same binning as used for the pythia Tune A MC. 
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FIG. 2: Measured detector level jet shapes in data for the tagged (grey full triangles) and inclusive (open squares) data samples 
for each pr bins. The jet algorithm used, MidPoint cone 0.7 j 'merge 0.75, is described in Sec. |III| 
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FIG. 5: Measured integrated 6-jet shapes for each of the four pr bins considered. The results are shown as black open squares 
where the error bars represent the statistical and total uncertainties. The statistical uncertainties are smaller than the squares. 
The results are compared to pythia Tune A (full lines) and herwig (dashed lines) predictions using the default /ib fractions 
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FIG. 6: This plot shows the same data points as in Fig. [5] The results are compared to pythia Tune A (full lines) and herwig 
(dashed lines) predictions if /ib is decreased by 0.2 (black lines). Also shown are the pythia Tune A and herwig predictions 
for single and double b-quark jets (light, highest, and dark, lowest, grey lines, respectively). 
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FIG. 7: Ratios of the hadron level integrated 6-jet shapes for various Monte Carlo predictions over the measured values. The 
light grey bands represent the total uncertainties on the measured b-jet shapes and the dark grey bands show the statistical 
uncertainties, pythia Tune A and herwig predictions using the default fib fractions are shown as grey lines (full lines for 
pythia Tune A, dashed ones for herwig). pythia Tune A and herwig predictions using f lh fractions 0.2 below the default 
values are also reported (black lines). 
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FIG. 8: Fractional pr outside a cone of radius r = 0.2 around the jet axis as a function of the pr of the jet. The results for 
6-jets are shown as black open squares and compared to different pythia Tune A and herwig predictions. The error bars on 
the plots represent the total and statistical uncertainties. The statistical uncertainties are smaller than the squares, (top) The 
results are compared to pythia Tune A (full lines) and herwig (dashed lines) predictions using the default /ib fractions (black 
lines). Also shown are the pythia Tune A predictions for the inclusive jet shapes (grey lines) as well as the previously published 
inclusive jet shape results (triangles), (bottom) The measured values are shown along with the expected distributions if /ib is 
decreased by 0.2 (black lines). Also shown are the pythia Tune A and herwig predictions for single and double fe-quark jets 
(light, lowest, and dark, highest, grey lines, respectively). 



